Relativistic Images in Randall- Sundrum II Braneworld Lensing 
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In this paper, we explore the properties of gravitational lensing by black holes in the Randall- 
Sundrum II braneworld. We use numerical techniques to calculate lensing observables using the 
Tidal Reissner-Nordstrom (TRN) and Garriga-Tanaka metrics to examine supermassive black holes 
and primordial black holes. We introduce a new way tp parameterize tidal charge in the TRN 
metric which results in a large increase in image magnifications for braneworld primordial black 
holes compared to their 4 dimensional analogues. Finally, we offer a mathematical analysis that 
allows us to analyze the validity of the logarithmic approximation of the bending angle for any 
static, spherically symmetric metric. We apply this to the TRN metric and show that it is valid for 
any amount of tidal charge. 
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I. INTRODUCTION 

Gravitational lensing provided the first experimental 
verification of General Relativity (GR) through obser- 
vations of starlight bending around the Sun during an 
eclipse in 1919 [H, i). Gravitational lensing is a very 
important probe of cosmological and astrophysical ques- 
tions [1,|4|. The vast majority of lensing studies are in the 
weak deflection limit, with light rays bent on the order of, 
at most, arcseconds. Over the last decade, there has been 
a renewed interest in strong deflection limit lensing and 
using it as a probe of GR. Some studies use a numerical 
method and a more analytical approach has been 

used as well In the last few decades, the mass of 

the supermassive black hole (SMBH) at the center of the 
galaxy (Sgr A*) has become well known [lj| E3, so the 
galactic center is often cited in the above lensing studies 
as an ideal candidate for the study of strong deflection 
limit lensing [H| QJ, El, El E3 ■ 

As the trajectory of a photon brings it closer a black 
hole's photon sphere, it undergoes a growing bending an- 
gle, and if the point of closest approach is incident with 
the photon sphere, this bending angle goes to infinity. 
This gives rise to a theoretically infinite sequence of im- 
ages close to the photon sphere on both sides of the optic 
axis due to photons looping around the black hole before 
reaching the observer. These images are termed relativis- 
tic images and their use as a probe of GR was pioneered 
by [H and expanded to include the topic of naked sin- 
gularities and cosmic censorship in 0, Q. However, as 
noted in these articles, relativistic images are highly de- 
magnified and their observational prospects are very dim 
at present. A later analysis [l2[ showed that a relatively 
"bright" relativistic image (brighter than the 32 nd mag- 
nitude in the K-band) may be visible in late 2017 through 
early 2018. Further investigations by [IHEillHl] discuss 
prospects for observation of secondary images that un- 



dergo a large bending angle. Studying the properties of 
these secondary images requires some of the techniques 
coming from studies in the strong deflection limit. As ob- 
servation of images lensed in the strong deflection limit 
is a distinct possibility, a close study of their properties 
is warranted. As mentioned above, 0, E2, HI] introduce 
a formalism that simplifes the calculation of observables 
for lensing in the strong deflection limit and introduces 
simple formulae for image positions and magnifications. 
This formalism can be used for any spherically symmetric 
metric and [l(J El, HH, H3 expand the analysis of strong 
field lensing to metrics which come from the Randall- 
Sundrum (RS) II braneworld scenario, a theory that has 
evoked great interest in the last decade. In Sec. HJ we re- 
view the two methodologies for solving the lens equation 
for lensing observables in the strong deflection limit, the 
numerical method of (H , and the analytical method based 
on approximating the bending angle with a logarithmic 
term by [Hill. In Sec. Hill we introduce the Randall- 
Sundrum braneworld and discuss several 4 dimensional 
solutions on the brane for the 5 dimensional black hole of 
RS II theory. We comment on several metrics that come 
from this theory and their regime of applicability. In Sec. 
TVl we calculate the observational properties of a galactic 
object lensed by the supermassive black hole at Sgr A*, 
and a galactic source lensed by a primordial black hole 
in our solar system. We compare the observables when 
modelling the black hole with the Schwarzschild metric 
as well as with several braneworld black hole metrics. In 
Sec. El we introduce a general test for the logarithmic 
approximation of the bending angle for all spherically 
symmetric metrics and show that it works for all values 
of tidal charge. Sec. |Vl] contains a discussion of the re- 
sults and possible future research directions. 



II. 



STRONG DEFLECTION LIMIT LENSING 
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In Sec. IJ we discussed the importance of lensing in the 
strong deflection limit. To go beyond the weak deflection 
limit for the bending angle, we will need to derive lensing 
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quantities directly from the metric. In this paper, we will 
only be considering metrics of the spherically symmetric 
form: 



ds 2 = -A(r)dt 2 + B[r)dr 2 + C(r)dn 2 



(1) 



The weak deflection limit approximation is only valid in 
the limit of a small bending angle and when the point 
of closest approach is far from the black hole. When the 
point of closest approach is close to the black hole, the 
bending angle is derived from the equations of motion 

M 



to be 



a(r ) 



( B{r) 
\C{r) 



1/2 



r a C(r) A(r ) 
V C(r ) A(r) 



dr 
r 



(2) 



-1/2 



(3) 



where r$ is the point of closest approach. Once we can 
calculate the bending angle of a null geodesic as a func- 
tion of ro, the results are used in conjunction with the 
Virbhadra-Ellis lens equation 



tan 8 = tan 8 — [tan ( 



+ tan(a - 9)] (4) 



where Dls, Dl, 0, and B are, respectively, the distance 
from the lens plane to the source plane, the distance from 
the observer to the lens plane, the angle of the image rela- 
tive to the optic axis, and the angle of the source relative 
to the optic axis. Ds, which is not pictured in Fig. [T]is 
the total distance from the observer plane to the source 
plane and is just Dls + Dl- In [23j . there is a discussion 
of lens equations used in gravitational lensing. While 
concluding that the "improved Ohanian" lens equation 
[H,[24[ is more accurate than Eq. ©, [23[ notes that the 
Virbhadra-Ellis lens equation is precise to within a factor 
of 1CP 4 in most situations. This suffices for the purposes 
of this paper and offers the best way to compare our re- 
sults with existent work, the majority of which is done 
with the Virbhadra-Ellis equation. While this is not an 
exact lensing solution [25j , this approximation is accurate 
and easy to use for calculating relativistic images [23| . A 
typical gravitational lensing scenario is pictured in Fig. 

m 

In a spherically symmetric scenario, the magnification 
of an image is given by 



sin/3 dB ._ 1 
( sm7 W 



(5) 



There are two distinct approaches in going from the 
bending angle to the observablcs of image positions and 




FIG. 1: The simplest system used to study gravitational lens- 
ing. A source at S, a lens at L, and an observer. 



image ma gnif ications. First we will review the approach 
of [ll|, I20I |26| that approximates the bending angle with 
a logarithmic term. The results in this section can be ap- 
plied to any spherically symmetric, static metric. Then 
we will review the completely numerical approach of 0- 
[§| in a Schwarzschild spacetime. This paper extends the 
numerical results in the literature to braneworld black 
holes. In a later section of the paper (Sec. fV| . we will 
show a general method for showing that the analytical 
approach reproduces the numerical approach for a par- 
ticular spacetime. 



A. Analytical Solution 

While an exact solution to the bending angle with- 
out reference to a background metric is given by pol \2l\ 
with solutions being given in integral form, the earliest 
attempt at trying to develop an analytic approach to rel- 
ativistic images was done by 0, where they attempted 
to express Eq. ((3]) as an expansion of an elliptic integral. 
They then found the first order expansion of the elliptic 
equation from its divergence at the photon sphere. This 
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approach is later used to calculate lensing observables in 
the braneworld scenario (Io| . The approach we will be 
working with in this paper was developed by [20I ] and we 
will develop the formalism with the intent of demonstrat- 
ing its range of effectiveness in Sec. |V] 

The equations of motions in Eq. ([2]) have cyclic coordi- 
nates t and <f> leading to conserved quantities 



E = A{r)i 
J = C{r)4> 

and the following expression for f : 



(6) 
(7) 



(8) 



where an overdot sepresents a derivative with respect to 
the affinc parameter. The radial and angular motion of 
the null geodesic can then be characterized by the ratio 
of the two conserved quantities: 




J 
E 



(9) 



In order for a photon which is initially travelling with 
r < 0, Eq. (|5j) must vanish for the photon to invert its 
motion and not fall into the black hole. For a given value 
of u and assuming the function ^ has a single minimum, 



will occur when 



C(r ) 



(10) 



where ro is the point of closest approach. Since § is 
lower-bounded, this equality can only be satsfied if u is 
greater than the minimum value 




where we have defined A m = A(r m ). The point r„ 
resents the radial coordinate of the photon sphere. 
These terms suggest the rewriting of Eq. ([3]) as 

B ,C 2\-1/2j 

C^A ~ ' 




(11) 



rep- 



(12) 



with the notation Di = Dls,E>l, breaking the bending 
angle integral into two parts- the first for the infall and 
the second for the portion of the photon's path from the 
radial minimum to the observer. To perform a detailed 
analysis of the bending angle around ro = r m , we will 
examine the function 



R(r,u) 



C{r) 
A(r) 



(13) 



From the previous discussion, we have shown that R{r, u) 
has a minimum at r m for any u. It also vanishes at 
(ro, it) by definition of tq and at (r m , u rn ) by definition of 



u m . We are interested in the properties of the bending 
angle corresponding to an inversion point very close to 
the photon sphere r rn , so we parameterize the inversion 
point as 



ro = r m (l + 8) 



(14) 



which corresponds to an impact parameter that is also 
very close to the minimum 



u m (l + e) 



(15) 



These properties allow for expansion of R{r, u) in orders 
of e and 8. The lowest order non-vanishing terms are: 



R(r,u) 



n 1 q2r , s 2 r2 9R, 

(16) 

and this creates the simple relationship between e and 8 
of 



Pm g2 



2ul 



where 



fin 



l d 2 R 
2F 



(17) 



(18) 



We have gone through this part of the derivation in 
detail because we will comment in Sec. |V] about the va- 
lidity of leaving off the higher-order terms in e and 8. 
From here, [20l | uses this expansion of R(r, u) to perform 
the integration in Eq. (Tj"2")) . Continuing along these lines, 
they define the position of the source {D^s) and the ob- 
server (Dl) with the coordinate r\ defined by: 



1-7? 



(19) 



with < 17 < 1. Eventually [20( finds that the bending 
angle is: 



A(j> = a\og 



S 2 



bo + bs 



(20) 



with the subscript O and S corresponding to the observer 
and source, respectively, and a and hi are functions of the 
metric 



En 



A m j3jj 



dr] Sign(rj)g{rj) 



(21) 
(22) 



g(v) = U r 



[R(T), U m ) 



-1/2. 




(23) 
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This result is equivalent to the result in [ljj when r\o = 
?/5 = 1. In order to calculate image positions and mangi- 
fications, we use the formalism of 



a 

a = — 
2 



b=-i 



g{r])dr) + alog^ 

-"m 



and adopt Eq. (9) from (Tlj 



n 

a{6) = Mog(^ - 1) 



(24) 
(25) 

(26) 



Using this logarithmic term in conjunction with a sim- 
plified lens equation 



P = B- ^Aa„ 



(27) 



where Aa n is the bending angle once 2im has been sub- 
tracted to account for the "loops" that the photon has 
made, we can easily solve for 9q which gives us the image 
position that corresponds to a bending angle a{9) = 2-nn: 



7T (1 + e n ) 



(28) 
(29) 



The image position for a given source position (3 is then 
given by 9q and a correction term: 



(0 - 0°) Dos 



(30) 



" " aD LS D L 

where the correction is much smaller than 8q. It is 
straightforward to derive the magnification of the n th rel- 
ativistic image as: 



u 2 m {l + e n )D os 
apU L L) LS 



(31) 



This outlines the analytical approximation in the strong- 
deflection limit. We will use it in this paper to briefly 
compare results with the purely numerical method. Also, 
we will examine if there are circumstances under which 
the higher order terms in Eq. (TTB")) would be relevant and 
therefore interfere with this form of the strong deflection 
limit approximation. We now review the numerical algo- 
rithm for gravitational lensing. 



B. Numerical Solution 

In the numerical procedure first outlined by [EH, con- 
servation of J and E is used to relate the image position 
with the point of closest approach. All the equations 
in this section up to here are true of any spherically 
symmetric, static metric. However, since numerical tech- 
niques will be employed in this section, the functions of 



the metric must be specified, and this section assumes 
the use of a Schwarzschild metric. In addition, this sec- 
tion and the rest of this paper uses geometric units in 
which c = G = 1 This approach can be applied to 
any metric, as it is in Sec. IIV1 The impact parameter 
u is conserved, and in the asymptotic Minkowski space 
(which is required for use of Eq. ©), the perpendicular 
distance from the center of the lens to the null geodesic 
is 



D L Sin 



(32) 



This is only valid if the source is in the asymptotic region 
(Dls >> 2M). At the point of closest approach, solving 
the geodesic equations yields 



2M\_i 
u = r (l ) 2 

TO 

To simplify, we use the following definitions: 



2M 
rp 

2M 



(33) 



(34) 
(35) 



to make the coordinates in terms of "Schwarzschild 
radii." 

Since Eq. (|32|) is equivalent to Eq. (|33]l. we can now 
describe the image position in terms of the coordinate of 
closest approach (as well as converting back from xq to 



Sin 9 



2M 



j_ 

J.'o 



(36) 



For a given value of /3 we insert Eq. ([3]) into Eq. d3| and 
solve for xq. This gives us a solution of the lens equation. 
Since Eq. (j4]) is transcedental, there are infinite positive 
and negative solutions corresponding to images on both 
sides of the lens. The solution with the highest (lowest) 
value of xo represents the classic lensing solution for the 
primary (secondary) image. To find a relativistic image, 
we must determine what range the bending angle would 
be in to yield a relativistic image. While it is possible 
that a source will align nearly perfectly with the optic 
axis and have a smaller angular position than the rela- 
tivistic image 0, in reality, a relativistic image will be 
closer to the optic axis (the line connecting the observer 
and source planes perp indicualar to both planes) than 
the source. In [l2|, [lj|, realistic possibilities for strong 
field images in the center of our galaxy do not occur with 
stars that are highly aligned with the optic axis (as the 
orbit of known stars near the center of the galaxy do not 
pass close to the optic axis). Hence, relativistic images 
will form at a smaller angle from the optic axis then the 
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a(xo) 




%0 Ana. 


A 


TT 


1.7603 


1.7084 


.02947 


2ir 


1.5451 


1.5433 


.00115 


3vr 


1.5091 


1.5090 


.00005 


47T 


1.5019 


1.5019 


2 x 10~ 6 


57T 


1.5004 


1.5004 


8 x 10" s 



TABLE I: The coordinate of closest approach required to yield 
a particular bending angle in both the numerical technique 
and the analytical approximation for a Schwarzschild space- 
time. The analytical approximation tends to be converge with 
the numerical technique as xo approaches the photon sphere, 
as the analytic approach loses validity as one moves away from 
the photon sphere. A is the fractional difference between the 
analytical and numerical result. See Sec. |V]for elaboration 
on the validity of the analytical approach. 

source angle and the bending angle for the first relativis- 
tic image on the side of the lens is slightly less than 2tt. 
Therefore, to find the correct bending angle that solves 
Eq. (01, the search for a correct value Xq should be in 
the region that yields a bending angle of 2tt. So, the first 
step is to compile a table of the values of xq that yield 
the values of mv. This is found in Table H To find a spe- 
cific relativistic image, be it higher order on either side 
of the lens, the bending angle must be estimated (it will 
be close to 2im for n th order images with sources close to 
the optic axis) and corresponding values for xq searched 
near the appropriate value. 

In addition to image position, image magnification is 
an important observational quantity. The general for- 
mula is given by Eq. ([5]), but incorporating Eq. (j4|) and 
expressed in terms of known functions of xq, magnifica- 
tion is given by: 



the small angle approximation can be utilized to make 
solving the equation more tractable, a point incorporated 
by Virbhadra into the numerical algorithm in [7|. One 
of the advantages of the numerical method is that it is 
universally valid and does not depend on approximations 
that are only valid in the strong or weak deflection limits, 
an important point when stud ying realistic scenarios for 
strong deflection limit lensing |l3l . Hij . 

III. BLACK HOLES IN THE BRANEWORLD 

An important and still unverified idea in modern 
physics is the existence of extra dimensions. This idea 
had its origins in the attempts of Kaluza and Klein to 
unify the electromagnetic force with GR by introducing 
an extra dimension, which was compactified in Klein's 
model (See reviews in [28], H^] ) . More recently, a class 
of these models was proposed where the Standard Model 
fields are constrained to a 3-brane while gravity prop- 
agates in a higher-dimensional bulk. Models which in- 
corporate ideas along this vein are colloquially referred 
to as brancworld models, starting with the outbreaking 
[30| . This paper will be based on a particular braneworld 
scenario of Randall and Sundrum [31 [32 1 . 

The Randall- Sundrum I model [3l| emerges from 11 
dimensional M-theory when the gauge fields of the stan- 
dard model are confined on two 1+9 branes located on 
the end points of an S 1 jZi orbifold, a scenario whose 
importance is highlighted by [33| and further developed 
by [3~4T - [3r| . 6 dimensions are compactified, making grav- 
ity effectively 5 dimensional. The Randall-Sundrum II 
scenario is this model with the second brane taken to in- 
finity [29j . The 5 dimensional metric in the absence of 
matter is known to be 1321: 



Sec' 



Dls 
D S 



[Sec 2 + Sec 2 (a(x o )-6)} 



Sec 2 /3 



da 



(37) 



Once we have solved for image position and Xq, the 
only unknown is 4§ which can be broken up into 



da da dxo 
dB dxo d6 



(38) 



and the calculation of these two quantities for the 
Schwarzschild metric is given in We have just il- 
lustrated how to go from a source position to an image 
position (given either as an angle, 6 or as a point of clos- 
est approach, Xq) and an image magnitude. 

One of the reasons why the numerical method is carried 
out is to make full use of Eq. (|4]) as the bending angle 
can be large and the weak field approximation does not 
hold. However, as Bozza points out, for a relativistic 
image, a is usually only slightly different than 27rn and 



ds 2 ^e- 2k ^[dt 2 -d 3 x]+dz 2 



(39) 



where k is the energy scale of the 5 dimensional cosmo- 
logical constant 



A, 



G 



-6k 2 



(40) 



and / is the characteristic size of the extra dimension. 
This cosmological constant prevents gravity from "leak- 
ing" into the bulk at low energies. The metric in the 
presence of matter is more complicated. Firstly, it is 
not clear whether a static black hole solution exists for 
a RS II braneworld. Although (37} has shown the exis- 
tence of a static black hole in a 2+1 brane setup, it is 
unclear what the solution is and whether there even is 
a static black hole solution in the corresponding 3 + 1 
brane scenario [38UI2I ] . One approach is to write down 
the induced Einstein equation on the brane and attempt 
to solve it. Following the derivation of [43|, the induced 
Einstein equations on the brane will be 

2 2 

= -Ag^ + k 2 T^ + - ^ + 4^- (41) 
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where A is the 4 dimensional cosmological constant which 
comes from a combination of the 5 dimensional cosmo- 
logical constant and the branc tension, k is related to the 
coupling constant in 5 dimensional as well as the brane 
tension A, S^ v is a term that is second-order in the stress- 
energy tensor, expresses contributions from the 5 
dimensional stress-energy tensor aside from the 5 dimen- 
sional cosmological constant [29|, |43| , and £ M „ is the dou- 
ble contraction of the Wcyl tensor with the unit normal 
to the brane. The derivation of Eq. (|41|) is performed in 
greater detail in (29j . To find the vacuum solution around 
a source, the following assumptions are made: The brane 
tension is finely tuned so that A in 4 dimensional is zero. 
Because we deal with a vacuum solution in 4 dimensional, 
the terms T^, J^"„„, and S, w are therefore zero outside 



/if! 



in' 



the source. This reduces the vacuum Einstein equations 
to 



(42) 



1, , . GM . 21 , 

-h 00 = V(r) «_(! + — 



(43) 



where /loo is a convention used for the difference from 
a Minkowski spacetime of the weak field metric term 
goo. Incorporating the weak field potential from Eq. (|43j) 
yields the following metric: 



dS\ = 



-(1- 



2M 4Mt 



2M 
-(1 + 

r 



3r 3 
2M1 2 
3r 3 



)dt 2 
)dr 2 



'■dn 2 



(44) 



This is the working metric for the weak field limit, 
or on scales r » I. Sometimes, these limits are not 
identical and this will be discussed in the next section. 
This metric is only expected to work in the linear regime, 
as the Ricci scalar does not vanish as it should (as ^ u is 
traceless) beyond first order in M [28| . 



In order to find a solution to the equations, some assump- 
tions must be made about the form of £ IIV . A weak-field, 
linear perturbative expansion is known |44(, but the so- 
lution has the wrong form in the strong field (28|. The 
non-linear strong field regime requires a different mctric- 
pcrhaps one metric will not be sufficient to cover the en- 
tire brane spacetime, but different metrics will apply in 
different regimes. Attempts to study the problem have 
yielded several possible black hole metrics [TsJ l44T - |47j and 
the effect of these metrics on lensing has been studied as 
well US EE llllli]- This paper will add to the literature 
by performing numerical studies of lensing using these 
metrics. 

The simplest solution, the "black string" was found by 
[46j which desribes an effective Schwarzschild solution on 
the brane (by setting the Weyl term £ M „ = 0). There arc 
a number of instabilities and difficulties with this solution 
[49| , and it is uninteresting from our perspective because 
the induced metric in 4 dimensions is Schwarzschild and 
observables in the Schwarzschild solution have already 
been calculated numerically. We will now set forth 3 
different metrics to be studied. 



A. Garriga-Tanaka 



B. Myers-Perry 

For small enough distance scales, we can consider the 
black hole as a five dimensional Schwarzschild black hole. 
On this length scale, the AdS curvature does not greatly 
affect the geometry of the black hole. Finding the metric 
can be done using using the Myers-Perry general form 
for the 4 dimensional induced metric for a higher dimen- 
sional black hole [5(1 [H| • This metric is [2l| 

dS 2 = -(1 - ^)dt 2 + (1 - ^r'dr 2 + r 2 dil 2 (45) 
where the black hole horizon radius rn is given by 




The approximate distance scale for which this metric 
can be considered accurate is r << I because that is 
the domain in which the graviton docs not differentiate 
between the AdS dimension and the other ones. lp and 
Mp are the known standard 4 dimensional Planck length 
and mass. 



An important solution for gravity in the braneworld 
scenario is the Garriga-Tanaka solution (44[. This met- 
ric comes from the fact that the linearized metric can be 
found at low energies (corresponding to scale r >> I). 
The perturbation on the brane is written in terms of a 
Green's function which is dominated by the low-energy 
zero mode of the 5 dimensional graviton. When calcu- 
lating the linearized correction to the metric for a point 
source, the Newtonian potential, to a first-order correc- 
tion is: 



Tidal Reissner-Nordstrom 



45] 



comes 



The "Tidal" Reissner-Nordstrom metric 
from the showing that all the contraints known for 
(symmetric, trace-free, and conservation equations) are 
satisfied by the stress energy tensor associated with the 
Einstein-Maxwell solutions of the standard 4 dimensional 
Einstein equations. By making the equivalence 

O (47) 



k 2 T„ 
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and using the appropriate choice of constants, a solution 
to Eq. dm is found: 



2M Q. ! 2 
(1 — — + ^) « r ' 



2 d0 2 
(48) 

Unlike the Reissner-Nordstrom solution, Q here is al- 
lowed to be negative, strengthening the effect of grav- 
ity. Since this metric is taken to apply only in the 
strong field, as its asymptotic behavior is incorrect for 
non-trivial values of Q, tidal charge is weakly limited by 
studies of neutron stars [52| and perhaps potential ob- 
servations of event horizon structure |53[ which are still 
inconclusive. The TRN metric is not favored for large 
black holes according to [4l[ , but it is an interesting can- 
didate for the strong field because the term encodes a 
5 dimensional potential that one would expect to see in a 
Randall-Sundrum braneworld. Studies [ljl[Il| about the 
lensing effects of the TRN metric around Sgr A* have set 
Q = q AM 2 where — 1 < q < 0, as this allows for signifi- 
cant effects from the braneworld term. 



IV. BLACK HOLE LENSING IN THE 
BRANEWORLD 

The greatest impetus for the study of relativistic im- 
ages, images formed deeped within the strong deflection 
regime near the photon sphere of a black hole, is the 
ability to test GR against an alternative theory of grav- 
ity. Photons that reach us after probing close to a black 
hole are most likely to carry signs of gravity's true na- 
ture. There have been previous studies of stron g de fl.ee- 
tion lensing around braneworld black holes [ToL 1 1 9L [2lj . 
However, all use some variant of the analytical formal- 
ism developed by [ll|, H3j ■ In this section, we will calcu- 
late lensing observables in two scenarios: The SMBH at 
the center of the galaxy lensing a galactic source and a 
small primordial black hole at solar system scales lensing 
a galactic source. For both, we will use the numerical 
technique. 

In order to choose an appropriate metric for calculat- 
ing relativistic image properties, we need to calculate the 
location of the photon sphere of a black hole. This is be- 
cause relativistic images need to have a very large bend- 
ing angle (at least tt for what [l(| terms "retro-lensing"). 
Such a large bending angle requires that the coordinate 
of closest approach for the null geodesic be very close to 
the photon sphere. Hence, knowledge of that scale tells 
us what regime we are dealing with and, therefore, which 
metric is appropriate. The photon sphere in braneworld 
metrics is given by setting (Eq. (11) in Q): 



2C(r)B(r)+r^^-B(r)-r^P-C(r) =0 (49) 
dr dr 



r% T = M 



3^ M 2 



(-5l 2 M + 3M 3 + V^VWlP+WW)^ 
(-5l 2 M + Ml 3 + V5V5Z 4 M 2 + 6/ 2 M 4 )^ 



3^ 



r MP 
ps 



V2r H 



(50) 
(51) 
(52) 



Note that these are all different than in the Scwarzschild 
metric in which r ps = 3M. We now follow the procedure 
outlined in Section II but use the appropriate braneworld 
metric. 



A. Supermassive Black Hole 

When considering the black hole at Sgr A* as a lens, 
we use the values D L = 8.3 kpc and M = 4.3 x 1O 6 M 
[l7j |. We model the source as a point source at D$ = 
2Dl. Modelling the supermassive black hole (SMBH) 
at Sgr A* as a braneworld black hole can be done with 
the Garriga-Tanaka metric, as the photon sphere is of 
order 10 7 km which is much greater than the scale of 
the extra dimension I. However, since the correction to 
the Schwarzschild metric is of order ~ 10~ 31 , the 
lensing corrections due to the Garriga-Tanaka metric is 
negligible as can be seen in Table [TTJ Although some 
numerical studies indicate that such a large black hole is 
unlikely to have an exterior metric significantly different 
than Schwarzschild [4l[ , a case can be made for the Tidal 
Reissner-Nordstrom metric as the spacetime curvature 
near the surface of the SMBH is much greater than on 
the surface of the Earth and, therefore, is not excluded by 
our tests of Newtonian gravity on Earth. The curvature 
of spacetime is measured by the Kretschmann scalar (as 
the Ricci scalar is zero in the vacuum outside a black 
hole), which is for the Schwarzschild spacetime 



K — RabcdR 



abed 



m 
Hr 



(53) 



The ratio of the Kretschmann scalar calculated at 
the surface (photon sphere) of the SMBH to the 
Kretschmann scalar calculated at the surface of the Earth 



KsgrA* 



2000 



(54) 



This yields 



This implies that gravity near the black hole is very 
strong and braneworld gravity may show up at in strong 
gravitational fields while not being detected on Earth. In 
addition, there are few constraints on positing a \ cor- 
rection term to the gravitational potential that only ex- 
ists in the strong field, so we can take the TRN metric as 
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a specific example of a ^ correction term to the metric. 
The TRN metric requires a choice as to the parametriza- 
tion of the variable Q. In the literature [lO, EH, Q 
is usually taken to have dimensions Q oc M 2 , so that 
the integral for the bending angle is dimensionless when 
rephrased in terms of Schwarzschild radii (See Eq. (|56|) ~). 
The dimensionless of the metric is desirable and therefore 
makes this parameterization of Q a natural choice: 



Q = q AM 2 



(55) 



with |g| usually taking on values between and 1. The 
weakness of this approach is that it makes the assump- 
tion that the strength of the bulk's "backreaction" onto 
the brane is oc M 2 . While this allows a neat result that 
eases calculations for relativistic observables and maxi- 
mizes the brancworld effect by maximizing Q in the case 
of large black holes, it is an ad hoc assumption. We will 
explore other paramcterizations of the backreaction onto 
the bulk for other cases of the black hole. For the case of 
the SMBH, we will use this parameterization as it gives 
the maximal amount of tidal charge. The results of the 
location of the Einstein rings for different values of q are 
contained in Table [TTJ Since the first relativistic image 
appears very close to the position of the first relativistic 
Einstein ring, and higher-order relativistic images appear 
close to their corresponding Einstein ring, the position 
and separation of the first two relativistic Einstein rings 
is a good indication of the properties of image positions 
for relativistic images. 

The effect of brancworld metric on image magnifica- 
tion can be seen in Tabic Mil Use of the TRN met- 
ric slightly demagnifies the images of a source almost 
directly directly behind (/? = 1.45/J,as) at the distance 
Ds = 2Dl- For the higher value of tidal charge, the 
images are even more demagnified. Repeating this cal- 
culation for multiple source positions ranging from 1 as 
to 1 /i as, the total magnification changes, but the ratio 
between the magnification in a Schwarzschild spacetime 
and a TRN spacetime for a specific value of q remains 
the same. For a SMBH, the GT metric will not produce 
any deviation from the standard GR results because the 
non- Schwarzschild term in the GT metric scales as ( — ) 2 . 
However, the TRN metric displays theoretically diffcr- 
cntiable results, a result that encourages further study. 
The images in the TRN results are fainter than the cor- 
responding images in the Schwarzschild metric for this 
particular case, but this is not at all a general charac- 
teristic of braneworld strong field images. This will be 
discussed further at the end of this section. 

Since we have chosen Q in Eq. (|48|) to be parameter- 
ized by q oc f the explicit form of Eq. ((3]) for the TRN 
metric when using the equality in Eq. (|35|) is: 



oitrn(xo) 



(It 



/jL(l_i + 4)_(l_ I + 4) 

V XQ V XQ Xq' \ x x A ' 

- it (56) 



Since there is no dependency on mass for the bending 
angle in either the TRN or Schwarzschild (q = 0) sce- 
narios, there is no dependency on mass in the ratio of 
the magnifications (Eq. ([57)) ) either. This is seen ex- 
plicitly by considering Eq. ([5|). As demonstrated by @, 
the position of relativistic images are very insensitive to 
source positions. Hence, the x coordinate position of the 
n th relativistic image is very close to the value of the x 
coordinate that yields a bending angle of 2im and the 
angular position of this image, 9, scales with M. If we 
hold (3 to be constant, the term ^ in Eq. (|5Tj) will scale 



as M. When we consider ^5, it is proportional to m dj5 
^ is constant in /3 (see graphs of image position against 
source position in [f|) and therefore the entire term scales 
as M. So the magnification of a relativistic image for a 
constant source position scales as M 2 . We can repeat the 
same argument for the bending angle in the braneworld 
scenario because Eq. (j56[) does not contain any depen- 
dency on M. Magnification would then have a scaling of 
M 2 as well. Hence, further studies of this form of the 
TRN metric will not yield interesting results, as the dif- 
ference between a Schwarzschild and a TRN lens does not 
change with the mass of the lens. In the next section, we 
will change the parameterization of the tidal charge for 
smaller black holes and anaylze the difference this makes 
for lensing properties. 



B. Primordial Black Holes 

One explanation for dark matter is primordial black 
holes (PBHs) that can form from a variety of mechanisms 
[Hi], [55| . Gould [56j j proposes femtolensing of gamma ray 
bursts, which are usually extra-galactical [601 ] . by com- 
pact dark matter. If PBHs make up a large fraction of 
dark matter in the universe, there is a chance of observing 
a gamma-ray burst with an interference pattern charac- 
teristic of lensing by a small black hole. An interference 
pattern resulting from such a lensing could potentially be 
observed by the ongoing Fermi Gamma-ray Space Tele- 
scope. This idea is extended to the braneworld black hole 
scenario by [48[ , who examine the effects of a GT metric 
of the interference pattern. A problem with analyzing 
a primordial black hole of the size considered by [48j |. 
1O~ 18 M0 is that relativistic images will pass in a very 
highly curved area near the black hole. The radius of 
curvature will be of order 10 -14 m, which is smaller than 
the wavelength of all electromagnetic radiation except for 
hard gamma rays, rendering the treatment unsuitable for 
the geometric optics approximation applied in this paper 
[3, ■ Therefore, we will consider a black hole of mass 
1CP 14 M0. In the Schwarzschild metric, the horizon size 
will be about 3 x 10~ u m and we can consider it to have 
some validity for even soft gamma rays. In the TRN 
spacetime, the horizon size will be even larger (about 
10~ 6 m for q = —.5), making our results applicable for 
the entire ultraviolet spectrum. Since the point of closest 
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Sch 


Garriga-Tanaka 


Tidal RN (q= -.1) 


Tidal RN (q= -.5) 


Einstein Ring 
1 st Relativistic Einstein Ring 
2 nd Relativistic Einstein Ring 


1.452 arcsec 
26.57 n arcsec 
26.54 /j, arcsec 


1.452 arcsec 
26.57 jj, arcsec 
26.54 jj, arcsec 


1.452 arcsec 
28.22 jj, arcsec 
28.19 /i arcsec 


1.452 arcsec 
33.38 jj, arcsec 
33.37 /i arcsec 



TABLE II: This table calculates the properties of Einstein rings formed by a source directly aligned with Sgr A* (M = 
4.3 x 1O 6 M0 and distance Dl = 8.3 kpc) and with a source distance of 2Dl- This table contains the angular positions of the 
primary Einstein ring and two relativistic Einstein rings formed by the source. The location of second relativistic Einstein ring 
is very close to the photon sphere and the change in its location in different metrics is a reflection of the different structure 
of the photon sphere in the braneworld scenario. Note the identical results for the Garriga-Tanaka and Schwarzschild metrics. 
Also note the difference in the gap between the first and second relativistic Einstein ring in each geometry. 





Sch 


TRN (q= -.1) 


TRN (q= -.5) 


Outer Rel. Image (opposite side of source) 
Inner Rel. Image (opp. side of source) 
Inner Rel. Image (side of source) 
Outer Rel. Image (side of source) 


-5.94 x 10" 12 
-1.10 x 10" 14 
1.10 x 10~ 14 
5.94 x 10~ 12 


-5.54 x 10" 12 
-7.85 x 10~ 15 
7.85 x 10" 15 
5.54 x 10~ 12 


-4.48 x 10" 12 
-3.97 x 10~ 15 
3.97 x 10" 16 
4.48 x 10~ 12 



TABLE III: This table contains the magnifications of relativistic images formed by a source at twice the distance of Sgr A* (see 
caption of Table [TTJ and at an angular position of /3 = 1.45 /xas or 10 -3 times the Einstein angle in the Schwarzschild geometry 
for a directly aligned source. These properties are calculated using several different metrics. The values for the GT metric are 
indistinguishable from the Schwarzschild metric and, therefore, have not been shown. 



approach is smaller than the scale of the extra dimension 
(r << I), I will ignore the Garriga-Tanaka metric and 
focus on the TRN and Myers-Perry metrics. 

While we can use the same parameterization as in Eq. 
([55]) . as explained before, the results would be qualita- 
tively similar. Also, it would be interesting to explore a 
scenario in which Q does not scale strictly as M 2 . In the 
paper introducing tidal charge, [45| considers a scenario 
in which tidal charge is fixed and does not depend on the 
mass of the gravitating body. For such a source of tidal 
charge, they cite an upper bound for Q of 



Q « 2M i? G 



(57) 



because higher values for Q will violate bounds on solar 
system tests if this metric is applied to our sun. This 
bound considers the possibility that Q is a fixed feature of 
geometry and does not scale with mass (this bound docs 
not apply if the TRN is proposed to only be the strong 
field limit of the metric). If we introduce this idea into 
our parameterization, the behavior of relativistic images 
in the braneworld scenario diverges from the behavior of 
Schwarzschild black holes. 

Now, I would like to consider a hybrid scenario. Let 
us say that we want Q cc M instead of M 2 , we would 
then need to introduce a mass scale in order to keep the 
metric and Eq. (|56p dimensionless. This will be a hybrid 
of tidal charge: part of it will come from mass getting 
reflected back at a linear rate and part will come from 
a fixed feature of the geometry. If Eq. (|57]l is to be 
satisfied, we must have M sca i e << Mq so I set it as 



then we can introduce the parameterization: 

Q = q 2M M scale 
This makes the bending angle integral: 



cxtrn(xo) 



(59) 



jLh- JL 

x V x 



qM 3calc 
IMxi 



qM a 



2M 



dx 



M, 



10~ 6 M ( 







Ira 



(58) 



It is evident from Eq. (JHUJ) that braneworld effects are 
inversely proportional to the lens mass. We will consider 
sources at solar system scales. We will model our lens 
as a primordial black hole at 1 AU and our source at 1 
Megaparsec. The source angular position /3 is set as one 
millionth of the Einstein angle for this configuration (in 
this case, 10 -12 as). In Table llVl we compare the dif- 
ference between relativistic images in the Schwarzschild 
and braneworld spacetimes. The large amount of tidal 
charge relative to the mass of the black hole magnifies 
images near braneworld black holes by several orders of 
magnitude compared to images near their Schwarzschild 
counterparts. 



C. Magnification in the Braneworld 

In the previous sections, we have shown that for the 
TRN metric, relativistic images in the braneworld are 
slightly fainter for relativistic images around supermas- 
sive black holes and are greatly enhanced for primordial 



(60) 
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Sch 


TRN (q= -.1) 


TRN (q= -.5) 


Myers-Perry 


Outer Rel. Image (opposite side of source) 
Inner Rel. Image (opp. side of source) 
Inner Rel. Image (side of source) 
Outer Rel. Image (side of source) 


-1.66 x 10 -29 
-3.08 x 10" 32 
3.08 x 10 -32 
1.66 x 10" 29 


-1.95 x 10 -21 
-2.69 x 10" 25 
2.69 x 10 -25 
1.95 x 10" 21 


-9.74 x 10 -21 
-1.35 x 10" 24 
1.35 x 10 -24 
9.74 x 10" 21 


-8.97 x 10 -25 
-1.24 x 10" 28 
1.24 x 10~ 28 
8.97 x 10~ 25 



TABLE IV: This table calculates the magnification of the outer two relativistic images on each side of the optic axis for a 
primordial black hole at 1 AU lensing a distant source at 1 Megaparsec, with /3 = 10 _12 /i as 



black holes. The distinction between the results in the 
case of primordial and supermassive black holes benefits 
from looking at the two components in Eq. ([5])- The two 
component magnifications are termed tangential magnifi- 
cation (/it) and radial magnification (/i r ) and are defined 



Sin 
H = (-, 



I'r 



Sm9' 
~ y d0> 



(61) 
(62) 



[ll| approximates both of these quantities and puts them 
in analytical form: 



A" 1 
V /in / 



aD L D LS ! 
Mr = (1 + 



(63) 
(64) 



From Eq. (|63|) . the tangential magnification of a rela- 
tivistic image is directly proportional to the image posi- 
tion of the relativistic image. This will always be larger 
when using a braneworld metric because of the strength- 
ened gravity in the braneworld metric. Since for primor- 
dial black holes, Q can be very large relative to the mass 
of the black hole, the photon sphere is approximated from 



Eq. ((521) as r, 



TRN 

ps 



when Q ^> M. So, the tan- 



gential magnification is always larger in the braneworld 
scenario, sometimes considerably greater. 

The variation in behavior of braneworld magnifications 
as mass changes comes from the radial magnification. 
From Eq. (|64p , it is apparent that for small radial magni- 
fications (which is always the case for relativistic images), 
radial magnifications are proportional to: 



fl r OC 



(65) 



This relation can be physically motivated- u m is a mea- 
sure of the photon sphere's size. This is because the coor- 
dinate of relativistic image is close to the location of the 
photon sphere. Magnification of the relativistic image is 
oc ^ , so when x is rescaled to reflect a larger photon 
sphere, the derivative becomes larger as well. The term 
e„ defines how quickly images decay towards towards the 



photon sphere for the n th image. If e n is greater, the im- 
age will shift more with a change in source position and 
therefore, the radial magnification will be brighter. The 
image brightness is also inversely proportional to a. 

It is simple to compute the behavior of each of these 
quantities in a braneworld scenario. In Fig. [5J we show 
the behavior of quantities a and b as q gets larger. The 
behavior of u m for large |<?| is given by Eq. (1ST " 



I -a 



-0.5 - 



-5 5 10 

Log (-q) 

FIG. 2: The behavior of strong field quantities as the behavior 
becomes 5 dimensional. The dotted line is the quantity a and 
the dashed quantity is b. From left to right, the quantity q 
becomes a progressively larger negative number. 



As can be seen from Fig. [2j the quantity a is equal 
to 1 for a Schwarzschild metric and becomes smaller as 
— q becomes large and asymtotes to The quantity 

b = —.400 in a Schwarzschild spacetime, but asymptotes 
to —.691 for a 5 dimensional spacetime. Therefore, e n 
(for n = 1) varies from .00125 in the Schwarzschild limit 
to .0000521 in the q — > — oo limit. When tidal charge is 
weaker (supermassive black hole) , u m is not much larger 
than the q = value, the smaller e„ term dominates and 
the relativistic image is demagnificd. When there is a 
great deal of tidal charge (Q ^> M), the great increase in 
u rn makes up for the e n term and the relativistic image 
is brighter compared to its Schwarzschild counterpart. 

For all black holes, relativistic images will be at a larger 
image position when using a braneworld metric, because 
the photon sphere of a braneworld metric is larger than 
its Schwarzschild counterpart. Hence, the tangential 
magnification is always higher using a braneworld. Ra- 
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dial magnification, as we have shown, depends on the 
value of the tidal charge. For a given source and lens 
position, the radial magnification is lower bounded, even 
as \q\ gets very large, but the tangential mangification 
is not upper bounded. For a large enough tidal charge, 
magnification will be enhanced for a braneworld image. 



V. APPLICABILITY OF THE ANALYTICAL 
METHOD 

In Sec. [TTl constructing a logarithmic approximation 
for the bending angle depends on the relationship in Eq. 
(|17[) between the impact parameter and the minimum ra- 
dial coordinate. When expanding R(ro,u) — in orders 
of e and <5, Eq. (|16[) is obtained by truncating at the first 
non-zero term in both r and u. However, the relationship 
between 8 and e is not being considered at an infinites- 
imal distance from the photon sphere. For example, in 
a Schwarzschild geometry, the first relativistic Einstein 
ring occurs with approximately ro ~ 3.09M Q, which 
corresponds to S = .03. This corresponds to e = —.0013. 
To be able to approximate the bending angle with a log- 
arithmic term, we must be able to neglect the higher 
order terms that do not appear in Eqs. (HHJ) and ([T7| . 
While [ill, HO] show that this first order treatment is ad- 
equate for typical spherically symmetric metrics such as 
the Schwarzchild, Reissner-Nordstrom (RN), and Janis- 
Newman-Winicour (JNW) metrics. However, some signs 
of this approximation breaking down is seen for RN and 
JNW metrics when they become near-extremal [ll| . Our 
goal in this section is to perform a general analysis of the 
higher order terms in the general case and then apply it 
to show the validity of the logarithmic approximation in 
the case of the TRN metric. To analyze the general case, 
we will expand out Eq. (|16p for one more order in both 
e and i5, which yields 



P m 5 2 + j m S 3 = 2w 2 „e + u 



2 £ 2 



where 



In 



ld 3 R 
6 dr 3 



(66) 



(67) 



First, consider the right side of Eq 



order 



term ^0-(r m ,u m ) is 
R(ro,u).We can neglect the e 2 term for the following 
reason. The ratio of the e terms is known- the ratio 
of the first-order term to the second order term is -. 

e 

This remains fixed regardless of the geometry because 
the u 2 term in R(ro,u) is independent of the metric. In 
the Schwarzschild case, e is much smaller than 6 and, 
therefore, higher order terms in e are less significant than 
higher-order terms in 5. This holds true in most space- 
times and can be checked using Eq. (|17|) . Examining the 
left side of Eq. (f66"]) . the ratio between the second order 
and third order terms in S is 



(Pl>. The higher 
from the definition of 



A, 



7m S 



(68) 



The greater this quantity is, the less significant the third 
order term is. Since /3 m and j m are functions of the met- 
ric, the ratio of thes two terms depends on the underly- 
ing spacetime. Since the validity of this approximation 
scheme is known for the Schwarzschild metric, it is useful 
to compare the ratio — in different spacetimes to its ra- 
tio in the Schwarzschild spacetime. For a Schwarzschild 
metric: 



ftn 
In 



(69) 



This result explains why the approximation scheme in 
Sec. |TT]only holds up for points of closest approach that 
are close to the photon sphere (small S). For large 5, 
the higher-order terms remain significant. However, as 
demonstrated in [ill . [20l | , it can be considered valid for 
the domain of relativistic images. Hence, we know that 
for the ratio in Eq. (jrJ§)) . we can consider the bending an- 
gle approximation to be valid. If, for an alternate space- 
time, this ratio becomes bigger, the approximation will 
be better. If the ratio becomes smaller, the approxima- 
tion will fare worse. For the TRN metric, the ratio — 
depends on the parameter q. As -q gets larger, the ratio 
gets smaller as well and the approximation gets worse as 
well. The analytical expression for both /3 m and j m can 
be obtained easily. Figure [3] displays the relationship be- 
tween q and ^ Fig. [3] shows that as q gets smaller, 




FIG. 3: The behavior of — as -q becomes larger. The de- 
crease in its value shows that the logarithmic approximation 
for the bending angle becomes worse as q gets smaller. 



drops from its starting value of | in a Schwarzschild met- 
ric. For the parameterization in Eq. (|59[) . the amount 
of tidal charge can be very large. We can evaluate the 
performance of the approximation in the limit of large 
tidal charge: 



,. Pm 1 

lim = — 

q^-co J m 3 



(70) 



This shows that for any amount of negative tidal 
charge, the approximation will be worse than in the 
Schwarzschild case, but only marginally so. This is borne 
out by an actual comparison of results obtained by both 
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methods. As S gets larger, the validity of the approxima- 
tion falls off quicker in the TRN spacetime. As q — > — oo, 
the TRN metric becomes a 4 dimensional metric, like 
the Myers-Perry metric, for which strong deflection limit 
lensing has been examined [13, HH . The result in Eq. 
(|70f shows that the strong deflection limit is accurate for 
the TRN metric and for the Myers-Perry metric. For any 
spherically symmetric metric, one test for this approxi- 
mation scheme can be the examination of — in that 

7™ 

spacetime. 

VI. CONCLUSION 

In this paper, we have shown various properties of rel- 
ativistic images and how the qualitative features of black 
hole lensing for a braneworld black hole depend on the 
size of the black hole. We have introduced a new pa- 
rameterization of "tidal charge" for the Tidal RN met- 
ric and have thereby uncovered a richer phenomenology 
of braneworld black hole structure in the strong deflec- 
tion limit. We show that the smaller the black hole is, 
the greater the effect of the braneworld scenario is on 
the relativistic image observables. The lensing proper- 
ties of small black holes are especially relevant with the 
start of the LHC, which has the possibility of producing 
TeV-scale black holes in the braneworld scenario. While 
the lensing properties of such a small black hole can- 
not be calculated using the geometric optics approach 
employed in this paper, there is good reason to suspect 
that lensing behavior in the braneworld will be dramat- 
ically different than Schwarzschild behavior. We have 
shown that relativistic images are greatly enhanced by 
the braneworld geometry, but it remains an open ques- 
tion whether there are observational consequences to the 
image enhancement. 

We also confirm that the analytical method serves as 
an accurate probe of relativistic images, within a reason- 
able margin of error. At this stage, when the observa- 
tional possibilitcs for relativistic images are more about 
the observation of the images than cataloguing the ex- 
act positions and magnitudes, a qualitative description 
of relativistic images is enough to explore this topic fur- 
ther. An advantage of the analytical approximation is 
that it allows for the calculation of certain observables 
more easily then using the numerical method. For exam- 



ple, the sum of the magnitude of all relativistic images 
after the first (which should appear unresolved) is eas- 
ily accomplished using the analyitical methodology. The 
analytical formulation also can be easily applied to ex- 
tended sources [l(|, which is more difficult with the nu- 
merical methodology. However, use of a numerical inte- 
gral is necessary in the astrophysically relevant situation 
p3 - tl4| of secondary images of S stars oribiting Sgr A* 
that undergo a bending angle in the middle of the region 
< a < ir. Here, both the weak and strong deflection 
approximations fail and evaluating integrals numerically 
is required for accurate calculation of primary and sec- 
ondary image properties. Hence, both methodologies will 
be relevant for the calculation of potentially observable 
consequences of relativistic images. 

It has been less than 10 years since the possibility of 
observing relativistic images and other strong field im- 
ages has been considered in and these images repre- 
sent a promising source for knowledge about gravity in 
its full strength. Observation of Sgr A* is a major area of 
research, and a theoretically interesting lensing observ- 
able can pave the way for attempts to use high-resolution 
VLBI arrays to find these elusive images. Upcoming in- 
struments such as GRAVITY and MICADO at the Euro- 
pean Extremely Large Telescope (E-ELT) also have the 
potential to make observations of strong deflection limit 
lensing near Sgr A*. Observation of a single relativistic 
image has the potential to answer deeply fundamental 
questions such as whether there is an extra dimension, 
scalar charge, or even black holes. The immense theoret- 
ical potential for relativistic observables in the vicinity 
of Sgr A* should motivate effort to continue to refine 
theoretical models and lensing scenarios. 

Acknowledgments 

I would like to thank R.K. Sheth for his unwavering 
support at all stages of research and writing. I thank 
K.S. Virbhadra for suggesting some aspects of this study, 
and help at the initial stages of this work. Thanks are 
due to J. Guzik, J. Hyde, J. Khoury and C. Dinerman 
for helpful discussions and to A. Congdon for reviewing 
the manuscript. I would also like to thank G. Cwilich 
and Yeshiva University for hosting me during research 
for this paper. 



[1] P. Schneider, J. Ehlers, and E. Falco, Gravitational 

Lenses ( Springer- Verlag, Berlin, 1992). 
[2] D. F. W., A. S. Eddington, and C. Davidson, Philos. 

Trans. R. Soc. London (1919). 
[3] J. Wambsganss, Living Reviews in Relativity 1 (1998). 
[4] H. Hoekstra and B. Jain, Annual Review of Nuclear and 

Particle Science 58, 99 (2008), 0805.0139. 
[5] K. S. Virbhad ra and G. F. R. Ellis, P hys. Rev. D 62, 

084003 (2000), [arXiv:astro-ph/9904193| 



[6] K. S. Virbhadra and G. F. R. Ellis, Phys. Rev. D 65, 
103004 (2002). 

[7] K. S. Virbhadra, Phys. Rev. D 79, 083004 (2009), 
0810.2109. 

[8] K. S. Virbhadra and C. R. Keeton, Phys. Rev. D 77, 
124014 (2008). 

[9] V. Bozza, S. Capozziello, G. Iovane, and G. Scarpetta, 
Gen. Rel. and Grav. 33, 1 (2001). 
[10] E. F. Eiroa, Phys. Rev. D 71, 083010 (2005). 



13 



V. Bozza, Phys. Rev. D 66, 103001 (2002). 
V. Bozza and L. Mancini, ApJ 611, 1045 (2004). 
V. Bozza and L. Mancini, ApJ 627, 790 (2005). 
V. Bozza and L. Mancini, ApJ 696, 435 (2009). 

E. F. Eiroa, G. E. Romero, and D. F. Torres, Phys. Rev. 
D 66, 024010 (2002), [arX"iv:gr-qc/0203049| 

F. Eisenhauer et al., ApJ 628, 246 (2005). 

S. Gillessen et al, astro-ph/0810.4674 (2008). 

V. Bozza, F. DeLuca, and G. Scarpetta, Phys. Rev. D 

74, 063001 (2006). 

R. Whisker, Phys Rev. D 71, 064004 (2005), 
|arXiv:astro-ph/041178"6| 

V. Bozza and G. Scarpetta, Phys. Rev. D 76, 083008 
(2007), 0705.0246. 

A. Majumdar and N. Mukherjee, International Journal 
of Modern Physics D (2005). 

S. Weinberg, Gravitation and Cosmology: Principles and 
Applications of the General Theory of Relativity (Wiley, 
New York, 1972). 

V. Bozza, Phys. Rev. D 78, 103005 (2008). 

H. Ohanian, Am. J. Phys. 55, 428 (1987). 

S. Frittelli and E. T. Newman, Phys. Rev. D 59, 124001 

(1999). 

V. Bozza, S. Capozziello, G. Iovane, and G. Scarpetto, 
Gen. Rel. Grav. 33, 1535 (2001). 

S. Frittelli, T. P. Kling, and E. T. Newman, Phys. Rev. 

D 61, 064021 (2000). 

R. Whisker, gr-qc/0810.1534 (2008). 

R. Maartens, Living Reviews in Relativity 7 (2004). 

N. Arkani-Hamed, S. Dimopoulos, and G. Dvali, Phys. 

Rev. Lett. B 429, 263 (1998). 

L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 
(1999). 

L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 
(1999). 

P. Horfava and E. Witten, Nuclear Physics B 460, 506 
(1996), |arXiv:hep-th/9510209| 

A. Lukas, B. A. Ovrut, K. S. Stelle, and D. Waldram, 



Phys. Rev. D 59, 086001 (1999), |arXiv:hep-th/9803235| 

A. Lukas, B. A. Ovrut, and D. Waldram, Phys. Rev. D 

57, 7529 (1998), [arX7v:hep-th/97lTT97l 

A. Lukas, B. A. Ovrut, and D. Waldram, Nuclear Physics 

B 532, 43 (1998), |arXiv:hep-th/9710208| 

R. Emparan, G. T. Horowitz, and R. G. Myers, Journal of 

High Energy Physics 1, 7 (2000), |arXiv:hep-th/9911043| 

P. Kanti and K. Tamvakis, Phys. Rev. D 65, 084010 



[40 
[41 

[42; 

[43 
[44 
[45 
[46 

I 47 

[48 

[49 

[50 

[51 

[52 

[53 
[54 

[55 

[56 
[57 

[58 
[59 
[60 



Classical and Quantum Gravity 23, 6633 (2006), 
|arXiv:hep-th/0606006| 

T. Shiromizu and M. Shibata, Phys. Rev. D 62, 127502 



(2000), [arXiv:hep-th/0007203 
T. Wiseman, Phy s. Rev. D 65, 124007 (2002), 
|arXiv:hep-th/0111057| 

A. 17. Fitzpatrick, L. Randall, and T. Wiseman, 
Journal of High En ergy Physics 11, 33 (2006), 
|arXiv:hep-th/0608208| 

T. Shiromizu, K. I. Maeda, and M. Sasak i, Phys. Rev. D 

62, 024012 (2000), [arXTv:gr-qc/9910076| 

J. Garriga and T. Tanaka, Phys. Rev. Lett. 84, 2778 

(2000). 

N. Dadhich, R. Maartens, P. Papadopoulos, and V. Reza- 
nia, Phys. Lett. B 487, 1 (2000). 

A. Chamblin, S. W. Hawking , and H. S. Reall, Phys. 

Rev. D 61, 065007 (2000), [arXTv:hep-th/9909205| 

R. Casadio, A. Fabbri, and L. Mazzacurati, Phys. Rev. 

D 65, 084040 (2002), [arXlv:gr-qc/011l"072| 

C.R. Keeton and A.O. Petters, Phys. Rev. D 73, 104032 

(2006). 

R. Gregory and R. Laflamme, Physical Review Letters 



(2002) , ||arXiv:hep-th/0110298| 

S. Creek, R. Gregory, P. Kanti, and B. Mistry, 



70, 2837 (1993), [arXiv:hep-th /9301052 

R. Guedens, D. Clancy, and A. R. Liddle, Phys. Rev. D 

66, 043513 (2002), |arX~iv:astro-ph/0205149| 

R. Guedens, D. Clancy, and A. R. Liddle, Phys. Rev. D 

66, 083509 (2002). 

A. Kotrlova, Z. Stuchlik, and G. Torok, Class. Quantum 

Grav. 25, 225016 (2008). 

S. Doelman et al., Nature 455 (2008). 

S. Hawking and B. Carr, Mon. Not. R. Astron. Soc. 168, 

399 (1974). 

D. Boyanovsky, H. de Vega, and D. J. Scwarz, Ann. Rev. 
Nucl. Part. Sci 2006. 56, 441 (2006). 
A. Gould, ApJ 386, L5 (1991). 

A. N. Aliev and P. Talazan, Phys. Rev. D 80, 044023 
(2009). 

F. De Paolis, A. Geralico, G. Ingrosso, and A. A. 
Nucita, Astron. A strophys. 409, 809 (2003), 
|arXiv:astro-ph/0307493] 

A. Petter, H. Levine, and J. Wambsganss, Singularity 
Theory and Gravitational Lensing (Birkhauser, Boston, 
2001). 

P. Podsiadlowski, P. A. Mazzali, K. Nomoto, D. Laz- 
zati, and E. Cappellaro, Astrophys. J. 607, L17 (2004), 
|arXiv:astro-ph/0403399] 



